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Abstract. This paper describes a Buchberger-style algorithm to compute a 
Grobner basis of a polynomial ideal, allowing for a selection strategy based on 
"signatures". We explain how three recent algorithms can be viewed as different 
strategies for the new algorithm, and how other selection strategies can be 
formulated. We describe a fourth as an example. We analyze the strategies 
both theoretically and empirically, leading to some surprising results. 



1. Introduction 

A fundamental tool of symbolic and algebraic computation is the method of 
Grobner bases. The first algorithm to compute a Grobner basis was introduced 
by Buchberger in 1965 [5]; subsequently, the computer algebra community has 
developed a number of additional algorithms, such as [4j [9j [T2l Q3] . In recent years, 
a new genus of algorithm has emerged, exemplified by F5 and G 2 V [TO] 111). 

While they are presented as different algorithms, both use a property called 
a "signature" to control the computation and reduction of 5-polynomials. While 
studying the two, we realized that they could be viewed as variations in the selec- 
tion strategy of a basic algorithm common to both. A third recent algorithm of 
Arri likewise fits this mold [2] , so it seemed instructive to formulate explicitly both 
the underlying structure and the three algorithms as strategies for implementation. 
This provides a common theoretical framework which allows a careful comparison, 
looking both at how the criteria employed by the strategies are related, and at 
experimental timings in a unified environment. For the latter, we employed both 
interpreted code (via Sage [15] and SINGULAR [7]) and compiled code (via Singu- 
lar). We examined both the original implementations and new implementations 
of the algorithms as "plugins" to the common algorithm. For consistency, the al- 
gorithms should all compute a reduced Grobner basis incrementally, so when the 
reader sees "F5", s/he should understand "F5C" [8]. 

Section [2] reviews basic notation and concepts, adapting different sources which 
vary considerably in notation (2J [8l QUI [Til OS]- Theoretical contributions, which 
include the new "sig-redundant criterion'', begin in Section 12.21 and continue into 
Section[3J Section 2] gets to the meat of comparing the strategies; both the analysis 
of Section 14.11 and the timings of Section 14.21 produce surprising and unexpected 
results. 

2. Background 

Let i G N, F a field, and R = F [xi, . . . , x n ]. Throughout this paper, F{ = 
• • • ) fi) where each fj e R. and Ii = (Fi) is the ideal of R generated by the 
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elements of F,. Fix a degree-compatible ordering < on the monoid M of monomials 
of xi, . . . , x n ; for any p G R, we denote p's leading monomial by lm (p), its leading 
coefficient by lc (p), and write It (p) = lc (p)lm (p). For brevity, we may denote 
t p = lm (p), c p = lc (p), and 1cm (t p , t q ) = t Viq . 

Let /i+i G R\Ii- We want an algorithm that, given a Grobner basis G; of Ii, 
computes a Grobner basis of Ii + i = (i^+i), where Fi + \ — (fi, . . . , fi+i). 

2.1. The traditional approach. Given p,q G i?, the S-polynomial of p and g 

is 

C Cp tpq 

Lp i^q Lq 

(It makes some things easier later if we multiply only q by a field element.) Further, 
given p,r G R and G C R, we say p reduces to r modulo G if there exist 
Ai, . . ., A.£ G N, h, . . . ,t£ G M, c%, . . . ,ce G F, and r , . . . , r e G R such that 

• r o — P an d ?'£ = r; and 

• for alH = 1, . . . ,£, 

• n = ri_i - CiUg At , and 

• lm (r ? ) < lm (n_i). 

Buchberger's algorithm computes the Grobner basis G of (-Pi+i) by computing 
S'-polynomials and reducing them modulo the current value of G: initially the 5- 
polynomials of /j, fi + \ where j = 1, . . . ,i; then, for any ^-polynomial that reduces 
to nonzero r, also for the pairs r, g where g G G, adding r to the basis after 
determining new pairs. Termination is guaranteed by Dickson's Lemma, since no 
polynomial is added to Gi+i unless it expands the M-submodule (lm (Gi+i)) of the 
Noetherian M-monomodule M |13| . 

The following is fundamental to the traditional approach. 

Definition 1. Let s G R and G C R, with #G = i. We say that s has a standard 
representation with respect to G if there exist hi , . . . , hi G R such that s = 
h\g\ + • • - + higi and for each k — !,...,£ either h^ = or lm (hk) lm (g^) < lm (s). 
We may also say that (hi, . . . , hi) is a standard representation of s with respect to 
G. 

If s reduces to r modulo G, then it has a standard representation modulo G; the 
converse, however, is often false. 

2.2. Signature-based strategies. All algorithms to compute a Grobner basis fol- 
low the basic blueprint of Buchberger's algorithm, but recent algorithms introduce 
a new point of view. We call these "signature-based strategies", inasmuch as their 
computations take "signatures" into account. 

Definition 2. Let Fi,...,F m be the canonical generators of the free R-module 
R m . Let p G (i*i+i), jeN with j <i + l, and hi, . . . ,hj G R such that hj ^ and 

P = hifi H hjfj. 

If c = lc (hj) and r = lm (hj), we say that crFj is a natural signature of p. Let 

§ be the set of all natural signatures: 

§ = {crFj : c G F\ {0} , r G M, j = 1, . . . , m} . 

We extend the ordering < on "ML to a partial ordering ~< on § in the following way: 
coFj -< drFfc iff j < k, or j = k and a < r. // j = k, a = r, and c = d, we say 
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that caFj = drFk- If j = k and a = r, we say that caFj and drFk are level. We 
do not otherwise compare them. 

Proposition 3. The ordering -< is a well-ordering on S = {tFj : r £ M, j = 1, . . . , m}. 

Thus, for eachp £ (Fi), uie can identify a unique, minimal, monic natural signature. 

(Our "Propositions" are either trivial or proved elsewhere.) 

Definition 4. We call the unique minimal monic natural signature of p £ R its 
minimal signature. We denote the set of all natural signatures of p by sig(p), 
and the minimal signature of p by S (p) . 

Proposition 5. Letp,q £ Assume that caFj £ sig(p) and drFk £ sig(g). 

Lett £ M. We have 

(A) caFj £ sig (p ± q) if caFj >- drF k ; 

(B) (c ± d) aF j £ sig (p ± q) if caFj and drFk are level and c ± d =/= 0; and 

(C) ctaFj £ sig(tp). 

Corollary 6. Let p,q £ ij+i. Assume that caFj £ sig (p) and drFk £ sig(g). 
Suppose that there exist a £ F and t £ M suc/i that at\t (q) = It (p). 

(A) If trFk -< o"Fj, f/ien ccrFj £ sig (p — atg). 

(^S^) // trFk, aF j are level and ad ^ c, then (c — ad) aF j £ sig (p — atq). 

Definition 7. If (caFj,p) , (drF fc ,a) £ § x I l+1 , a £ F, and t £ M satisfy (A) or 

(B) of Corollary® we say that p — atq is a a-reduction of p with respect to 
q. Otherwise, p — atq is a-unsafe. When it is clear from context that we mean 
a-reduction for appropriate a, we refer simply to reduction. 

We define a a-reduction of p modulo G analogously, and say that it is com- 
plete when no reductions of the type described in Corollary are possible. It is 
semi-complete when reductions of type (B) can be performed, but not reductions 
of type (A). 

We now adapt the notion of a standard representation to consider natural sig- 
natures. 

Definition 8. Let G C Sx/, +1 with #G = I, and suppose that for each (drFj,g) £ 
G we have tF j = S (g). We say that any (ccrF,--|_i, s) £ § x I i+ i has a standard 
representation with respect to G (or sig-representation, or a -representation 

for short) if caFi + \ £ sig (s) and there exist hi, . . . ,h# £ R such that (hi, . . . , hi) 
is a standard representation of s and each hk = or crFi + i >^ lm (hh) S (gfc). 

Just as a reduction of p to zero modulo G corresponds to a lm (p)-representation 
of p with respect to G in the traditional case, a cr-reduction to zero modulo G 
corresponds to a cr-representation with respect to G. 

Theorem 9. Let c\, . . . , q £ F, ri, . . . ,rt £ M, and pi, . . . ,pt £ {0} such 



G = {{a 1 F jl ,g 1 ),..., (a M F jM , g M )} 

U {(F i+ i,/ i+ i) , (diTiF i+ i,pi) , . . . , (dtT £ F l+1 ,p e )} , 

aF j = S (g) for all (caFj,g) £ G, and for all (co"Fj_|_i,p), (drFj,q) £ G one of the 
following holds: 



that 
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— • a \ Fi + i >- ( -f 3 - ■ t) Fj and ( -j 3 - • (jFj_|_i, S Pt q J /ias a standard represen- 



ts 

tation with respect to G. 
Then G = {g : 3 (o~Fj,g) £ G} is a Grobner basis of 

Proof. Recall from [3] that G is a Grobner basis of 7,-_|_ i iff ^G^) = /j+i and S^g 

has a standard representation with respect to G for all distinct p,q 6 G. Since 
Gi is a Grobner basis of h, we know that S p , q has a standard representation for 
every p, q £ Gi, so it suffices to check only the pairs p,q £ G such that at least 
p G;. For any such pair where (-f 3 - ■ cm F,*_|_i >- ^-jr 2 - ■ rj Fj, by hypothesis 

f "f* 2, • crFi-)_i, S'p.q^j has a standard representation with respect to G; by definition, 

5 Pl9 also has a standard representation with respect to G. 

Order the remaining ^-polynomials by ascending level signature, and choose 
one such S'-polynomial S p , q such that (caFi + i,p) and (drFi + i,q) are in G and 

(if 3 - ■ a) F i+ i = (tz* ■ tJ F i+ i is minimally level. Now, s = ^ ■ p - % • ^ ■ q 

has signature less than (jf 3 - • o~\ F i+1 ; let hi, . . . ,hj £ R such that s = \Y^ k=1 hkfk 

and \m{hj)Fj = S(s). Each lm (hk) lm (/ft) Fft is smaller than (jf 3 - ■ cr^j F i+ i, 
the mimimal level signature. By hypothesis, iS-polynomials corresponding to top- 
cancellations among these lm (hk) lm (/&) have sig-representations with respect to 
G. Thus, there exist Hi, ... , -H#g such that s = J] HkQk, and for each fc, iJ^ = 
or lm (Hk) lm (gk) < lm (s) and 5 (s) >z hn (i/fc) 5 (gk)- If s = Sp jg , we are 
done. Otherwise, lm (s) = t v . q , implying lm (He) lm (ge) = t VA for some I. Thus, 
there exist u £ M and (fj,Fj,g) £ G such that ubxi(g) = t Ptq , and (uu)Fj -< 

(jf 3 - ■ crj F i+ i. Since (uu)Fj -< (jj 3 - ■ crj F i+ i, by hypothesis S^g and 5 9l9 have 

sig-representations with respect to G; since S^.^ = -f^Sp.g + j^Sg^q, so does S Pl q. 

For the remaining 5-polynomials, proceed similarly by ascending level signature. 
Top-cancellations of level signature will be smaller than the working signature, and 
so will have been considered already. □ 

We need one more concept. It looks innocent, but is quite powerful. We have not 
seen it elsewhere in the literature, but it is inspired by other work; see Lemma [T"5l 

Definition 10. Let G C SxR and (ccrF i+1) /) £ SxR. If there exists (drF i+ i,g) £ 
G such that r | a and lm (g) | lm (/), then (crFj+i,/) is signature-redundant to 
G, or sig-redundant for short. 

The proof of Theorem [9] uses a technique common to signature-based algorithms: 
proceed from smaller to larger signature, reusing previous work to rewrite S'-poly- 
nomials. This suggests an algorithm to compute a Grobner basis; see Algorithm Q] 
It follows the basic outline of Buchberger's algorithm, with several exceptions. 

• As we will show in Lemma II2[ Algorithm [T] prepends a natural signature 
to each critical "pair", and considers pairs by ascending natural signature 
(as in mm 03] [16]), rather than by ascending 1cm (as in [HH]). 

• Only cr-reductions of (<7Fj_|_i,r) are computed directly. We require semi- 
complete reduction, but complete reduction implies this. If we perform a 
complete reduction and conclude with (ccrFj+i, r), we multiply r by c _1 to 
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Algorithm 1 Signature-based Grobner basis computation 



l: inputs 

2: Fi C R, such that Fi is a Grobner basis of (Fi) and fj ^ • • • , fj— l) 

3: £ i?\ (Ji) 

4: Outputs 

5: G C § x R satisfying the hypothesis of Theorem [9J 

6: do 

7: Let G = {(F u fx) , . . . , (Fi, h) , (F, +1 , / <+1 )) 

8: LetP={(^fF i+1 ,f i+ x,f):feF i } 

9: Initialize 5yz - TBD 

10: while P ^ do 

11: Prune P using %2 TBD 

12: Let S — {(o"Fj-|_i,p, q) 6 P : deg cr minimal} 

13: Let P = P\S 

14: while S ^ do 

15: Prune S using Syz, G — TBD 

16: Let (aFi+x , P, o) G <5 such that crFi + i is minimal 

17: Remove (crF.j + i,p, q) from £ 

18: Let (cFj-i-i, r) be a semi-complete cr-reduction of S Pl q mod G 

19: Update fusing (aF i+1 ,r) -TBD 

20: if r and (crFj+i, r) not sig-redundant to G then 

21: for (rFi + x,g) <E G such that g ^ and g not sig-redundant do 

22: if • (7 ^ • r then 

23: if • cr > ■ t then 

t r tg 

24: Let (/xF i+ i,j5,g) = (jff- ■ o-F i+1 ,r,s) 

25: else 

26: Let (>F i+ i,p,eO = (^-rF^^.rj 

27: if deg [i~d then 

28: Add (/iF <+ i,jj, g) to S 

29: else 

30: Add (/j,F i+ x,p, q) to P 

31: Append (crFi + i,r) to G 

32: return {((jF, + i,j)£G: j/0, not sig-redundant} 



ensure crFi + i £ sig (r) for line [THJ Reductions that are cr-unsafe occur in 
line 1251 via the generation of new critical pairs. Algorithm [1] adds these to 
S rather than P to preserve the strategy of ascending signature. 

• The if statement of line [50] rejects not only zero polynomials, but sig- 
redundant polynomials as well. This has a double effect in Lemma Q2] and 
Theorem Q3J 

• We adopted the following from F5, partly to illuminate the relationship with 
this algorithm better. Algorithm [1] is easily reformulated without them, in 
which case it begins to resemble G 2 V and Arri's algorithm. 

• Critical pairs are oriented: any (crF; + i,p, q) £ P corresponds to S p , q = 
up — cvq where (rF,_|_i,p), (fJ>Fj, q) € G and cr = ur >- u/x. 
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• Line [12] selects all pairs of minimal degree of natural signature. With 
homogeneous polynomials and a degree-compatible ordering, this se- 
lects all .^-polynomials of minimal degree. An inner loop processes 
these by ascending signature. 

• Zero and sig-redundant polynomials are retained in the basis for rea- 
sons that become clear later; however, line [21] prevents them from 
being used to compute new critical pairs, and line [32] does not add 
them to the output. 

Remark 11. For now, we define lines [P[ [77] 1 15\ and \19\ to do nothing, and discuss 
them in Section \2.S\ 

We prove the correctness of Algorithm [JJ in several steps. 

Lemma 12. Suppose that one of lines [R \28l or\3(A creates (crF j +1 , p, q) . Write 
s = S p „; n °t on ly is £ s ig( s )t but ^i+i = S (s) unless s already has a 

sig-representation w.r.t. G when Algorithm^ would generate it. 

Proof. That crFi+i £ sig (s) follows from Proposition [5] Corollary[6j and inspection 
of the algorithm. For the second assertion, suppose that crF;+i ^ S (s). Let rFj = 
S (s); by definition, there exist hi, ... , hj such that j < i + 1, s = hif± + ■ ■ ■ + hjfj, 
hj 7^ , and lm (hj) = r. Notice rFj -< crF, + i. Since the algorithm proceeds 
by ascending natural signature, top-cancellations of smaller signature have been 
considered already. Hence, all top-cancellations among the h^fk would have sig- 
representations at the moment Algorithm [1] would generate s. We can therefore 
rewrite the top-cancellations repeatedly until we conclude with a sig-representation 
of s. □ 

Lemma 13. Suppose that line \20A prevents the algorithm from creating critical pairs 
using (crF.; + i,r). Then r = 0, the corresponding S -polynomials have sig-represen- 
tations already, or will after consideration of pairs queued in P U S . 

Proof. If r = 0, then we are done. Suppose r ^ 0; by line[20l there exist (rFj+i , g) € 
G such that r | a and lm (g) \ lm(r). Let u € M such that ulm (g) = lm(r). 
If ut < a, then line [TS] did not perform a semi-complete er-reduction of S p _ q , a 
contradiction. Hence ut > a. 

Let t £ M such that ut > a = ir, so lm (r) = ulm (<j) > tlm (g). The sig- 
nature /iFj of s = r — tg is smaller than ctF.; + i, so Algorithm Q] has considered 
top-cancellations of this and smaller natural signature. Hence, (fj,Fj,s) has a [i- 
representation. Critical pairs have been generated for g, so for any (£F; + i, q) g G, 
S r ^ q = ur — cvq = u (s + tg) — cvq, whose top-cancellations already have a sig-repre- 
sentation or will after consideration of pairs queued in P U S. □ 

Theorem 14. Algorithm^ terminates correctly. 

Proof. Correctness: If we show that the output of the algorithm satisfies the hy- 
pothesis of Theorem then we are done. By Lemma [T^] and the strategy of 
ascending signature, we know for any (crFj,g) £ G that g = or aFj = S (g). 
The only S'-polynomials for which the algorithm does not explicitly compute sig- 
representations are those satisfying the criteria of the if statement of line [20] and 
the criterion of line [22] The criterion of line 1201 is the hypothesis of Lemma [TBI with 
it and the criterion of line [22] we complete the hypothesis of Theorem [9] 
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Termination: Let M' be the monoid of monomials in xi, . . . , X2n', as with M, we 
can consider it to be a Noetherian M'-monomodule. Any (o"F,+i,r) added to G 
with r 7^ corresponds to an element of M' via the bijection 

(<7,im(r))= (n^%n^) ii V, 

Let J be the M'-submodule generated by these elements of G. Suppose the algo- 
rithm adds (crFj-|_i,r) to G and J does not expand; this implies that there exists 
(rFj+i, g) E G such that t \ a and lm (g) | lm (r). Since (o\EV)_i, r) is sig-redundant, 
line 1201 prevents it from generating new pairs. 

Hence, every time Algorithm [T] adds (crF, + i,r) to G, either the submodule J 
expands, or the algorithm abstains from computing pairs. A submodule of M' can 
expand only finitely many times, so the algorithm can compute only finitely many 
pairs. Hence, the algorithm terminates. □ 

The following interesting result will prove useful; its criterion is used in [TT| 
to prevent the generation of new pairs. 

Lemma 15. To see «/ (oF-j+i, r) is sig-redundant in Algorithm^ it suffices to check 
if there exist (rF^+i, g) £ G and t £ M such that tr = a and ttm (g) = lm (r). 

Proof. Assume that there exists (rFj+i, g) <E G such that r | a and lm (g) | lm (r). 
Let t, u £ M such that tr = a and itlm (g) = lm (r). If ut < a, then line IT8l did 
not compute a semi-complete cr-reduction of S Piq , a contradiction. If ut > a, then 
u > t, so lm (r) = ulm (<?) > tlm (g). The signature of r — tg is smaller than er, so 
r — tg has a sig-representation with respect to G. In addition, lm (r — tg) = lm (r); 
by the definition of a sig-representation, there exist (fjIFj, h) £ G and u € M such 
that ulm (h) = lm (r — tg) = lm (r) and u ■ /iFj is no greater than the signature of 
r — tg; that is, u ■ [iFj -< <tFj_|-i. But then line [18] did not compute a semi-complete 
cr-reduction of S Piq , a contradiction. □ 

2.3. Pruning P and S. This section does not propose any criteria that have not 
appeared elsewhere; rather, it lays the groundwork for showing how lines I§1 HT1 [T5l 
and [TH] can use such criteria to improve the efficiency of Algorithm [1] The general 
idea is: 

• Syz will consist of a list of monomials corresponding to known syzygies; i.e., 
if t G Syz, then £F.; + i is a natural signature of a known syzygy. 

• Line [TT1 removes (crFj+i,p, q) from P if there exists t £ Syz such that t | a. 

• Line [15] does the same, and ensures that if (crFj-|_i,p, q), (oFj+i,/, g) € S, 
then at most one of these is retained. 

Already, Lemma [T2l suggests: 

(NM) Discard any (oFj + i, /, g) G P U 5 if oFi + i is not the minimal signature of 
Proposition 16. Line\Q can put Syz= {lm (<?) : <? £ Fj}. 

For a proof, see Lemma 16 in [8] (Faugere's Criterion). It is similar to the proof of 
the following criterion [5] [TT] : 

Lemma 17. If the result of UneUM is (oFj+i, r) r = 0, £/ien line{T^ can add 
a to Syz. 
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(Sketch). Suppose that line [18] gives (oF,-|_i,r) with r = 0. Now, r is the a- 
reduction of s = S p ^ q from line !18l and by Lemma [12] crFi+i £ sig (s). By definition, 
3h\, . . . , /li+i € R such that s = Y>hkfk and lm (/ij+i) = er. Since r = 0, there exist 
Hi, . . . , H#g such that s = ^H^gh, each fffc = or lm (-ff/c) lm (gk) < lm(s), 
and S (YlHkgk) ~< crFj+i- Hence T,hkfk — 51Hk9k = and has natural signature 
aFj_|_i. Suppose there exist (rFj_|_i,p, (?) G P U 5 and u £ M such that uer = r; 
then S'p.g = Sp.q — u(T,hkfk — ^Hkgk) has signature smaller than r; now apply 
Lemma [T^] □ 

Another criterion is implied by the following lemma. 

Lemma 18. Let r G M, B — {(ojFi+i, /j) G G : <jj \ r}. We can choose any 
((tF,-|_i,/) G i? and discard in line \15\ any (rFj_|_i, p, g) E P if p =/= f , or if we 
compute (tF.;_|_i, /, g) where p = f and g =/= q. 

(sketch). Choose any (crF^+i, /) G B, and let (/iFi+i,p) , (/i'F,;+i, q) G G such that 
(rFj+i,p, g) £ PUS. Let t,u G M such that ta = ufj, = t. The signature of 
up — t/ is smaller than r, so it has a standard representation with respect to G; 
say up - tf = hkgk- Then = up - vq = u (tf + h k gk) - vq. All top- 
cancellations in this representation of S Pjq are of equal or smaller natural signature, 
so S Pt q will have a standard representation with respect to G once line [TBI chooses 
o-Fi+i >- rF i+1 . □ 

Notice that Lemma [TBI requires only divisibility; if there are no (rF,; + i,p, q) £ S 
such that p = /, then we could discard all (rFi_|_i,p, q) G S. We thus have a 
"rewritable" criterion: 

(RW) For any rFi + i G § select (oFj_|_i,/) G G such that cr | r; discard any 
(rFj_|_i , p, q) if p ^ /, or if p = / and we retain another (rFi_|_i , /, g) G S 
where q ^ g. 

For simplicity's sake, we assume that we apply (RW) only in line [15] but (NM) 
both there and in line [TTJ 

3. Known strategies 

This section sketches briefly how Arri's algorithm, G 2 V, and F5 can be viewed 
as strategies for Algorithm [JJ distinguished by: 

(1) whether reduction is complete or semi-complete; and 

(2) how they prune P and 5*. 

Space restrictions prevent us from going too far into each algorithm's workings, 
or proving in detail the characterization of each as a strategy for Algorithm [JJ 
However, the reader can verify this by inspecting the relevant papers. 

3.1. Arri's algorithm. Algorithm [JJ is very close to Arri's algorithm, which uses 
semi-complete reduction. Although [3] presents this algorithm in non-incremental 
fashion, with a more general way to choose the signatures, we consider it incremen- 
tally, with the definition of signature as given here. 

The algorithm maintains a list G similar to that of Algorithm [JJ and discards 
Sf_ g if /, g do not satisfy a definition of a "normal pair". This differs from the 
definition in |10] : 

Definition 19. Any f,g G are a normal pair if Sf. g = uf — cvg and 
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• for any (uFj+i , p) G {(S (f) , /) , («S (g) ,g)} there does not exist (rFj+i, q) £ 
G and ( g M such that tr — a and tlm (q) = lm (p); 

• S (it/) = u ■ S (/) and S (vg) = v ■ S (g); and 



In addition to G, Arri's algorithm maintains a list L of leading monomials used 
to prune P (there called B). These correspond to known syzygies; whenever s 
sig-reduces to zero, the monomial part of its natural signature is added to L. 

We can characterize this as a strategy for Algorithm [1] in the following way. The 
first bullet of Definition [T9l implies the sig- redundant property (Lemma [T5)l . To im- 
plement the second bullet, [2] counsels initializing L to {lm (/) : / G F,} and adding 
a to L if the cr-reduction of r concludes with 0. This implements Proposition [16] 
and Lemma 1171 In addition, [2] points out that for any fixed natural signature 
one should keep a polynomial of minimal leading monomial; after all, cr-reduction 
occurs when the leading monomial decreases and the natural signature is preserved. 
Thus, the algorithm discards any S'-polynomial if another polynomial of the same 
natural signature has lower leading monomial. This implements Lemma 1181 So, 
Arri's algorithm discards (crFi+x,p, q) if either of the following holds: 
(AM) for some g G F i: lm (g) | cr, or for some (rF i+1 , r) G G, r | a and r = 0; or 
(AR) there exist (rF i+1 , g) G G and t G M such that tr = a and lm (tg) < 
lm (S Pt q), or there exist (rFj+x, /, g) G S U P and ( 6 M such that tr = a 
and lm {tSf y3 ) < lm (S p>q ). 
Notice that (AR) checks divisibility of a, not equality. 

Proposition 20. Arri's algorithm implements Algorithm]^ with semi-complete re- 
duction: (AM) implements (NM) and (AR) implements (RW). 

3.2. G 2 V. Although G 2 V can be used to compute the colon ideal, we consider it 
only in the context of computing a Grobner basis. Thus, we are really looking at a 
special case of G 2 V. 

G 2 V maintains two lists of polynomials, U and V. The polynomials of V are the 
elements of the basis. The polynomials of U are paired with those of V such that 

• if Vj G Fi, then Uj = 0; 

• if Vj = fi+i, then Uj = 1; and 

• if vj = ctkVk ± dtgvg for some c, d G F, tk,ti G M, and Vk,ve G V, then 
Uj = ctkUk ± dt£U£. 

In addition, 5-polynomials and reductions are computed in such a way that lm (iij) 
is invariant for all j: G 2 V computes Vj —ctvk only if lm (uj) > t\m (uk) or lm (iij) = 
tlm (iik) but lc (uj) lc (cuk)- Thus, if Uj ^ 0, then UjFi + i G sig (vj). 

The algorithm maintains another list H of monomials that is initialized with the 
leading monomials of all f £ Fi, and expanded during the course of the algorithm by 
adding lm (iij) whenever Vj reduces to zero. It does not compute an S'-polynomial 
for Vj and Vk if: 



In addition, if every possible reduction of Vj is by some Vk such that It (vj) = 
dt\t (vk) and It (uj) = dvlt (uk), then Vj is super top-reducible, and [11] abstains 
from generating critical pairs for S Pl q if: 

(GS) either p or q is super top-reducible. 



• S(uf)^S(vg). 
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Criterion (GM) implements Proposition [16] and Lemma |T7l while (GS) implies 
the sig-redundant property (Lemma IT5j) . G 2 V offers no implementation of (RW) 
beyond, "store only one [pair] for each distinct [natural signature]". Which pair is 
left somewhat ambiguous, but we will see that the choice is important. 

Proposition 21. G 2 V implements Algorithm [JJ with complete reduction: (GM) 
implements (NM). 

3.3. F5. As explained in the introduction, we use the F5C variant of F5 [5]. In 
fact, we actually use a simplified version of F5; we describe the differences below. 

F5 maintains several lists G±, . . . , Gj+i C S x R; for j = 1, . . . , i, each Gj is a 
Grobner basis of (/i, . . . , fj). Whenever a (<tF,+i, r) concludes er-reduction, (a, r) 
is added to the (i + l)-st list in a list named Rules. 

F5 discards (crFj + i,p, q) if: 
(FM) for some g S Fj, lm (g) | er, or 

(FR) there exists (t, g) £ Rules i+ i, not sig-redundant, such that c/ was computed 

after p and r \ a . 
Notice that (FR), like (AR), checks divisibility of a, not equality. 

Proposition 22. The simplified F5 described here implements Algorithm [JJ with 
semi-complete reduction: (FM) implements (NM), and (FR) implements (RW). 

As noted, the F5 described here is simpler than |10| . where: 

• S'-polynomials of minimal degree are not computed in any particular order 
(but the code of j8] proceeds by ascending 1cm rather than ascending natural 
signature) ; 

• if S p _ q = up — cvq, then Criteria (FM) and (FR) are applied not only to up 
but to vq; in addition, for any potential cr-reduction r — tg, the criteria are 
used to reject some tg. 

Omitting these does not represent a significant difference from the original algo- 
rithm. In fact, descriptions of F5 by ascending signature have been around for some 
time f[TI I16|): the second bullet can be viewed an optimization that makes sense in 
an F4-style implementation, such as [TJ. 

There is one significant difference: the original F5 would not check for sig- 
redundant polynomials. In view of this, when we view F5 as a strategy for Al- 
gorithm [TJ we will include the sig-redundant criterion; but in Section 14.21 we will 
look at F5 with and without this criterion. 

That said, in Section 1431 we will look at both the simpler F5 described here, and 
an implementation of the original. 

4. Comparison of the algorithms 

The thrust of Section[3]was to show that Arri's algorithm, G 2 V, and a simplified 
F5 can be viewed as implementations of Algorithm [TJ Section 14.11 by contrast, 
compares the three algorithms as their authors originally defined them, using a 
strictly logical comparison of which critical pairs are discarded, without regard to 
timings. Nevertheless, we retain the notation of Algorithm [TJ 

Section 14.21 compares the three algorithms both ways: as implementations to, 
or "plugins" for, Algorithm [JJ and for G 2 V and F5 as standalone implementations. 
The results between the approaches (plugin vs. original) do not differ significantly, 
but provide both surprising results and additional insights. 
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4.1. Logical comparison of the algorithms. In this section we consider care- 
fully how the criteria of Arri's algorithm, G 2 V, and the simplified F5 overlap. We 
begin with complete vs. semi-complete reductions: 

Fact 23. G 2 V reduces using some polynomials that F5 and Arri's algorithm do 
not. 

Proof. This is because reductions are complete in G 2 V, but only semi-complete in 
F5 and Arri's algorithm. □ 

On the other hand: 

Fact 24. In Algorithm^ there cannot exist a reduction of type (B) in Corollary^ 
without a reduction of type (A). 

Proof. Let (rFj+i, r) , (crFj+i, g) E G such that lc (g) = d, lc (r) = c ^ d, and 
there exists t E M such that lm (r) = t\m (g) and r = to. Then r — tg has a 
natural signature smaller than rFj+i; since the algorithm proceeds by ascending 
natural signature, r — tg has a sig- representation. Since lc (g) ^ lc (r), we have 
lm (r — tg) = lm (r), so the sig- representation of r — tg implies that there exist 
u E M, (/xF, + i, /) E G such that vim (/) = lm (r) and ufi < t. This is a reduction 
of type (A). □ 

Together, Facts [^51 and [^1 mean that while G 2 V can reduce a polynomial using 
polynomials that F5 and Arri's algorithm cannot, they can still reduce it using 
other polynomials. 

The next observation regards non-trivial syzygies. 

Fact 25. (GM) and (AM) are equivalent. In addition, some pairs (cPi+i , p, q) 
rejected by (GM) and (AM), but not by (FM), are also rejected (FR). 

Proof. It is trivial that (GM) and (AM) are equivalent. Inspection shows that (FM) 
=> (GM), but not the converse. Assume therefore that (GM) rejects a critical pair 
(crFj_|_i,p, q) that (FM) does not; this implies that a is divisible by some r E M, 
where (tF,+i , r) was the result of a complete or semi-complete reduction, and r = 0. 

Fact [2U implies that r tr-reduces to zero in F5 as well. This is recorded in 
Rules by appending (r, r) to Rulesi + i. Since r = 0, F5 generates no more critical 
pairs for it. So p, q ^= r. If r was generated after p, then Rewritten(cr) =/= p, so 
Rewritten?(crFi + i,p) would return True. In this case, (FR) rejects (rFj+i,p, g). 
Hence at least some pairs (erF i +1 , p, q) rejected by (GM) but not by (FM) are also 
rejected by (FR). 

On the other hand, suppose that F5 computed p after r; then Rewritten(crFi + i) =/= 
r. In fact, it might return p, and Rewritten?(erFi + i,p) would return False, so 
that F5 would not reject (<rF i+ i,p, q), whereas (GM) would. □ 

On the other hand, Lemma [T7l implies that one could modify (FM) to consider 
zero reductions as well as trivial syzygies, as the other algorithms do. We try this in 
the next section. Alternately, one could modify (FR) to scan Rulesi + \ for pointers 
to zero reductions, using them to discard pairs before performing the usual (FR) 
criterion. Either works, but the former would likely be more efficient in interpreted 
code; see a related discussion in the following section. 

Fact 26. Some critical pairs computed by F5 are not computed by G^V and Arri's 
algorithm. 
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Test case 


F5 


G 2 V 


Arri 


Arri+(MR) 


Katsura-9 


14.98 


17.63 


18.25 


20.95 


Katsura-10 


153.35 


192.20 


185.76 


220.01 


Eco-8 


2.24* 


0.49 


0.45 


0.53 


Eco-9 


77.13* 


13.15 


5.20 


14.59 


Schrans-Troost 


3.7 


5.3 


6.46 


7.72 


F744 


19.35* 


26.86 


7.37 


27.77 


Cyclic-7 


7.00 


33.85 


8.82 


40.54 


Cyclic-8 


7310 


26242 


17672 


>8h 



*See the discussion in the text. 
Table 1. Timings, in seconds, of compiled Singular implemen- 
tations of the strategies, implemented as plugins for Algorithm [TJ 
Computed on a workstation with a 2.66GHz Intel Core 2 Duo 
P8800 and 4 GB RAM, running 64-bit Ubuntu 10.10. Base field is 
IF32003; ordering is degree reverse lexicographic. 

Proof. From Lemmall5[ we know that G 2 V and Arri's algorithm do not compute 
critical pairs for sig-redundant polynomials, whereas F5 does. □ 

Fact 27. Some pairs (erF , p, q) discarded by (FR) are not discarded by G 2 V. 
Likewise, some pairs (erFj_|_i,]9, q) discarded by (AR) are not discarded by G 2 V. 

Proof. As noted in Section l3~2| G 2 V implements (RW) by checking for equal signa- 
tures only, whereas (FR) and (AR) check for divisibility. As a consequence, (FR) 
and (AR) can discard ((7Fj_i_i, p, q) because (rFj+i,/) G G and r | cr, even if / 
generates no pairs of natural signature a. □ 

4.2. Experimental results. Although Facts [25] and [26] imply an advantage for 
G 2 V and Arri's algorithm over F5, Fact [27] implies an advantage for F5 and Arri's 
algorithm over G 2 V. We will see that the latter advantage is more significant than 
the former. 

We first look at some timings of the algorithms as plugins for Algorithm [1] To 
do this, we implemented Algorithm [I] as a C++ class in the kernel of a devel- 
oper version of Singular 3-1-2, then created descendant classes corresponding to 
the other algorithms. This allowed us to implement complete reductions for G 2 V 
and semi-complete reductions for F5 and Arri's algorithm without giving either an 
otherwise unfair advantage. Using compiled code allows us to avoid the overhead 
of an interpreter, but the code was otherwise unoptimized, linking to Singular's 
polynomial arithmetic. Table [T] lists timings in seconds corresponding to this im- 
plementation. We do not include timings for a bare-bones Algorithm [TJ having no 
criteria to prune P or S, it is unbearably slow. In Table [2] we count the number of 
critical pairs reduced, along with the number of zero reductions. 

As an example, we considered a fourth strategy that is essentially Arri's algo- 
rithm, but we replace (AR) by 

(MR) there exist (rFj + i ,g) G G and t G M such that tr = a and g has fewer mono- 
mials than S Pt q, in which case we consider tg in place of S p ^ q ; if (cr, f,g)(zS 
then we may choose either (cr,p, q) or (cr, /, g) freely. 
This implementation of (RW) causes more work than necessary. The first polyno- 
mials generated tend to have the fewest monomials, so (MR) selects these instead 
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Test case 


F5 


G^V 


Arri 


Katsura-9 


886;0 


886;0 


886;0 


Katsura-10 


1781;0 


1781;0 


1781;0 


Eco-8 


830;322* 


2012:57 


694:57 


Eco-9 


2087;929* 


5794;120 


1852;120 


Schrans-Troost 


380;0 


451;0 


370;0 


F744 


1324;342* 


2145;169 


1282;169 


Cyclic-7 


1063;44 


3108:36 


781:36 


Cyclic-8 


7066;244 


24600;244 


5320:244 



*See the discussion in the text. 
Table 2. Number of critical pairs reduced by each strategy, fol- 
lowed by number of zero reductions. We omit Arri+(MR) for 
space; it is comparable to G 2 V. 



of later polynomials, and so repeats many earlier reductions. The algorithm still 
computes a Grobner basis, but takes the scenic route. 

In general, F5 terminated the most quickly, but there were exceptions where 
Arri's algorithm did. This is explained by the discussion in the proof of Fact [25] 
(FR) is sometimes too aggressive, and does not notice some zero reductions. We 
modified (FM) to check for these first, and this modified F5 terminates for Eco-8 
(-9) in 0.38s (8.19s), and for F744 in 8.79s. Regarding critical pairs, it computes 
565 (1278) critical pairs for Eco-8 (-9), of which 57 (120) reduce to zero; and 1151 
critical pairs for F744, of which 169 reduce to zero. This suggests that computing 
lm (S Pt q) in order to check (AR) is usually too expensive for the benefit that we 
would expect, but in some cases it may be worthwhile. 

The results of Table Q] surprised us, in that it contradicts the unequivocal as- 
sertion of QT] that G 2 V is "two to ten times faster" than F5. Apparently, this 
is because compared implementations and not algorithms. Why is this prob- 
lematic? Primarily it is due to the use in [11] of interpreted code. Some natural 
adjustments to the implementation of F5 used in [8] [TT| are needed merely to start 
making the two comparable. We tried the following: 

• Formerly, the F5 implementation checked (FM) using an interpreted for 
loop, but the implementation of G 2 V checked (GM) using Singular's 
reduce (), pushing the for loop into compiled code (line 173, for exam- 
ple). We changed the F5 code to check (FM) using reduce (). 

• As specified in |10j . TopReduction and CritPair return after each reduc- 
tion or creation of a critical pair. This back-and-forth incurs a penalty; it 
is sensible to loop within these functions, returning only when reduction is 
semi-complete or all critical pairs have been considered. (The code for G 2 V 
did this already.) 

• In the code accompanying (8], Spol computes S-polynomials by ascending 
lcm, but one could proceed by ascending natural signature instead [1] I16| . 
(The code for G 2 V did this already.) 

Some other changes contributed a little; see Table [3] for timings, and Section [5] for 
source code. With this new implementation, G 2 V sometimes outperforms F5, but 
by a much smaller ratio than before. Tellingly, F5 outperforms G 2 V handily for 
Cyclic-n. This is despite the persistence of at least one major disadvantage: the 
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F5/G Z V 
in EH 


Test case 


F5 


G"V 


F5/G"V 


Katsura-6 


1.2 


0.53 


2.26 


6.32 


Katsura-7 


12.2 


6.7 


1.82 


4.91 


Katsura-8 


134.28 


50.1 


2.68 


5.95 


Schrans-Troost 


261.69 


50.61 


5.17 


14.04 


F633 


16.19 


2.72 


5.95 


14.50 


Cyclic-6 


7.1 


7.3 


0.97 


3.90 


Cyclic-7 


693.3 


962.5 


0.72 


3.12 



Table 3. Timings, in seconds, of Singular libraries, not imple- 
mented as pfugins to Algorithm [1] Computed on a MacBook Pro 
with a 2.4GHz Intel Core 2 Duo and 4 GB RAM, running OS X 
10.6.5. Base field is F32003; ordering is degree reverse lexicographic. 



code to check (FR) still uses an interpreted for loop. This imposes a penalty not 
only in Spol but also in IsReducible (called f ind_reductor in the implementa- 
tion); the interpreted for loop checking (FR) is one reason IsReducible consumes 
about half the time required to compute a Grobner basis in some systems. 

For G 2 V, on the other hand, we observed a disadvantage inherent to the al- 
gorithm and not to the implementation: its implementation of (RW) checks only 
equality, not divisibility. Fact [27] is especially evident in Table [2] Most 5-poly- 
nomials are subsequently discarded as super top-reducible. These discards are not 
reflected in Table 2 of |11| . which mistakenly implies G 2 V computes fewer polyno- 
mials than F5. Time lost in reduction is unrecoverable. 

A final note. When we first implemented G 2 V we obtained much worse timings 
than those of Table [TJ When we inspected the source code accompanying [11] , we 
found that the choice of which pair to store for a given signature is not arbitrary: the 
implementation prefers the most recently computed polynomial that can generate 
a given signature (lines 176-183). That is, it discards (<jFj + i , p, q) G S if 

(GR') there exists (crFj+i, f,g) € P such that p ^ / and / was computed after p, 
or V — f an d q was computed after g. 

Note the similarity with (FR). 

5. Conclusion 

This paper has described a common algorithm for which F5, G 2 V, and Arri's 
algorithm can be considered strategies, implemented as "plugins". Algorithm [1] 
makes use of a new criterion to reject polynomials and to guarantee termination, 
not only for itself, but for G 2 V and Arri's algorithm as well, through Lemma [T5l 
The matter is not so clear for F5, which seems to terminate all the same; we have 
not yet determined if some mechanism in F5 implies the sig-redundant criterion. 

Both timings and logical analysis imply that claims in [TT] that G 2 V is "two to ten 
times faster than F5" are based on a flawed comparison of implementations, rather 
than criteria. Indeed, the inherent advantages appear to lie with F5 and Arri's 
algorithm: in cases where F5 is not the most efficient, Arri's is, and modifying 
(FM) to check for non-trivial syzygies usually turns the scales back in F5's favor: 
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(FM') max y-f 3 - ■ a, ^j 3 - ■ t J is divisible by lm (g) for some g £ Fi, or by fi for some 

OuFi+i,0) eG." 

Source code developed for Table [3l along with a demonstration version of the algo- 
rithms for the Sage computer algebra system |15) . can be found at 

www . math . usm . edu/perry/Research/ 

by appending to the above path one of the filenames 

f 5_ex . lib, f 5_library . lib, f 5_library_new. lib, or basic_sigbased_gb .py . 
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